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A group theoretic approach to shear-free radiating stars 
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A systematic analysis of the junction condition, relating the radial pressure with the heat flow in 
a shear-free relativistic radiating star, is undertaken. This is a highly nonlinear partial differential 
equation in general. We obtain the Lie point symmetries that leave the boundary condition invariant. 
Using a linear combination of the symmetries, we transform the junction condition into ordinary 
differential equations. We present several new exact solutions to the junction condition. In each 
case we can identify the exact solution with a Lie point generator. Some of the solutions obtained 
satisfy the linear barotropic equation of state. As a special case we regain conformally flat models 
which were found previously. Our analysis highlights the interplay between Lie algebras, nonlinear 
differential equations and application to relativistic astrophysics. 


I. INTRODUCTION 

Radiating stars in a general relativistic context have 
been widely studied because of their importance in as¬ 
trophysics. The detailed physics of such models have 
been investigated by Herrera et al. |H with dissipation, 
by Di Prisco et al. [1] with charge, and Herrera et al. 
[3] with viscous dissipation in casual thermodynamics in 
the streaming and diffusion approximations. Particular 
models have been recently studied by Chan et al. !5[ for 
an imperfect nonadiabatic distribution, by Pinhciro and 
Chan I-}; with charge, and Pinheiro and Chan Q for a col¬ 
lapsing body with an initial inhomogeneous density. Her¬ 
rera and Santos 0 introduced the concept of Euclidean 
stars in general relativity by requiring that the areal ra¬ 
dius and proper radius are equal. Exact models satisfying 
this condition were generated by Govender et al. [§] and 
Govinder and Govender :9j. Abebe et al. 0 showed 
that generalized Euclidean stars may be generated with 
an equation of state. An interesting feature of their ap¬ 
proach is the application of Lie symmetries to produce 
models which are expanding, accelerating and shearing. 
This indicates that the Lie symmetry approach is helpful 
in solving the boundary junction condition, and high¬ 
lights its role in astrophysical applications. We point out 
that in recent times Lie symmetry generators have been 
particularly useful in generating models in spherically 
symmetric gravitational fields. Msomi et al. [Tl], 1121 have 
considered applications with heat flow and generated new 
classes of exact solutions. Kweyama et al. 00 anal¬ 
ysed Noether and Lie symmetries for charged relativistic 
fluids and were able to solve the field equations exactly. 

The shear-free assumption in stellar models is of¬ 
ten made in the study of relativistic radiating stars. 
Kolassis et al. 0 obtained the first exact solution 
with dissipation effects, satisfying the boundary condi- 
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tions, for a shear-free radiating star in geodesic motion. 
Thirukkanesh and Maharaj 0] generated exact solutions 
for the geodesic model by transforming the junction con¬ 
dition into Bernoulli, Riccati, and confluent hypergeo¬ 
metric equations. This model was extended by Govender 
and Thirukkanesh 0 to include a nonzero cosmological 
constant. Note that in these treatments fluid particles 
travel along geodesics and the Friedmann dust solution 
in the absence of heat flow was regained in all cases. 
An important point here is that these treatments were 
performed under the condition of pressure isotropy and 
anisotropic pressures are absent. Ivanov 0] considered 
a general shear-free perfect fluid with heat flow that con¬ 
tains conformal flatness and geodesic models as a special 
case. A conformally flat relativistic model with dissipa¬ 
tion and inhomogeneity was first proposed Herrera et al. 
0. This was integrated by Maharaj and Govender |20j 
and Herrera et al. M to obtained classes of exact so¬ 
lutions. Another class of conformally flat solutions was 
generated by Misthry et al. 0] by transforming the 
junction condition equation to an Abel equation. A fur¬ 
ther conformal flat radiating model, generated by a Lie 
symmetry was obtained by Abebe et al. 0] . Other mod¬ 
els of shear-free radiating stars have been constructed by 
Tewari 0 , 0 , Pant et al. 0 and Pant and Tewari 
0 which describe massive radiating fluid spheres and 
generated horizon-free collapse. 

The shear-free assumption is utilised often in the anal¬ 
ysis of self-gravitating spheres and when investigating 
gravitational collapse. This is largely due to the resulting 
simplification of the field equations. In addition, shear- 
free and the homogeneous expansion rate conditions are 
equivalent to the homology conditions in the Newtonian 
limit. Therefore the shear-free condition is well justi¬ 
fied. However we need to make the observation that the 
shear-free condition may be unstable in the presence of 
anisotropic pressures and dissipative fluxes. The condi¬ 
tions under which an initially shear-free fluid continues 
to remain shear-free in its subsequent evolution has been 
studied by Herrera et al. 0. In that analysis it was 
demonstrated that pressure anisotropy and dissipation 
affect the propagation of time and the gravitating rela- 
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tivistic fluid can become unstable. 

We analyse a shear-free radiating model, without as¬ 
suming conformal flatness and geodesic motion of fluid 
particles, in the presence of pressure anisotropy. We con¬ 
sider the Einstein field equations and boundary condi¬ 
tions in the context of Lie symmetries with the objective 
of producing new models of relativistic radiating objects. 
This is the main object of this paper. The case of shear¬ 
ing stellar models with Lie symmetries has been consid¬ 
ered in some earlier treatments. Those treatments do 
not allow for vanishing shear and our results cannot be 
regained from such investigations. The absence of shear 
produces fundamentally a different set of field equations 
to be integrated. The Lie generators are generically dif¬ 
ferent when the shear vanishes. In Sect. [©we briefly dis¬ 
cuss the shear-free radiating star and present the junction 
conditions. We show that the master equation is a non¬ 
linear partial differential equation in the metric functions. 
In Sect. m we generate Lie symmetries for the master 
equation. Using the Lie point symmetries approach we 
transform the boundary condition to ordinary differential 
equations. By analyzing the resulting ordinary differen¬ 
tial equations and transforming in the original variables 
we present new exact solutions in Sect. HVHVl andlVIl We 
link our solutions to specific Lie symmetry generators. In 
Sect. IVIII we show that some of our solutions satisfy a 
linear barotropic equation of state. We make concluding 
remarks, discuss the causal heat transport equation and 
summarize our solutions in tabular form in Sect. EEED 


II. THE MODEL 


We consider the particular case of spherically symmet¬ 
ric, shear-free radiating stellar models. The line element 
for the interior for the spacetimes is given by 

ds 2 = -A 2 dt 2 + B 2 [dr 2 + r 2 ( dd 2 + sin 2 ddcj) 2 )} , (1) 

where A and B are metric functions of t and r. The 
acceleration and and expansion are nonzero but the fluid 
is shear-free. 

The energy momentum tensor has the form 


Tab = ( a * + P±) u a u b + p±g ab + ( P\\ - P±)XaXb 

+q a u b + q b u a , (2) 


with heat flux and anisotropic stress. The fluid four- 
velocity u a = -^6 g is comoving, y a is an unit four-vector 
along the radial direction (u a x a = 0), and the heat flow 
vector q a = (0, c/,0,0) is radially directed ( u a q a = 0). 
The Einstein field equations for the heat conducting 


spherically symmetric anisotropic fluid © become 
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for the line element ©■ The equations ([3|) describe the 
gravitational interactions in the interior of a shear-free 
spherically symmetric star with heat flux and anisotropic 
pressures. 

The boundary of a radiating star divides the spacetime 
into interior and exterior regions. The interior spacetime 
© has to match across the boundary of the star to the 
Vaidya spacetime 

ds 2 = _ dy 2 _ 2dydR 

+R 2 ( dd 2 + sin 2 9d(/) 2 ) , (4) 

which is the exterior. Here the quantity m{v) denotes the 
mass of the star as measured by an observer at infinity. 
Matching leads to the junction conditions 
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where the hypersurface E defines the boundary of the 
radiating sphere. The particular junction condition (I5dl) 
should be solved, together with the field equations ©, to 
obtain the potentials A and B. This completes the model 
of a relativistic radiating star. The junction condition is 
a nonlinear differential equation 
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valid at the boundary of shear-free radiating star. Equa¬ 
tion © is the master equation that governs the evolution 
of the model. 
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III. THE MASTER EQUATION 


Using the invariants © we can write © in the form 


We can use the Lie analysis of differential equations to 
find a solution to the master equation ©. In a general 
relativistic context, the Lie analysis has proved useful in 
cosmological settings as shown by Msorni et al. fill . [T3 , 
Kweyama et al. [131 [lij and Nyonyi et al. [29l l30( |. and 
in astrophysical applications as illustrated in the treat¬ 
ments by Abebe et al. EI3, mi. The general approach 
is described in these treatments and we will not repeat 
details here. Essentially we have to find infinitesimal Lie 
point generators that allow us to reduce a partial differ¬ 
ential equation to ordinary differential equations. This 
is made possible by the existence of invariants associated 
with the Lie point symmetries of the partial differential 
equation equation ©. 

The process has been adapted over time and is now al¬ 
gorithmic, and so can be implemented by different com¬ 
puter software packages. Utilizing the program PRO¬ 
GRAM LIE 32], we find that © admits the following 
Lie point symmetries: 


Gi 

G 2 

g 3 


-Af'(t) 

A dA 
d 


A 


d_ 

dA + 
d_ 
dB 
d 


m 


dt ’ 


B 


dA dr ’ 


(7a) 

(7b) 

(7c) 


[2 a 2 gx 2b+2c+1 ((6 + c)g — xg')] h'h 2 
+2acx b+c+2 g 2 g' h' h — 2 c 2 x 2 g 3 g'h' 

+ [c 2 xg 2 (xg' 2 - 2 g((b + c- 1 )g' - xg"))] h 
- [ 2 acgx b+c+1 (g ( xg" + g') - xg' 2 )] h 2 
+ a 2 x 2< - b+c ^ (xg 1 — bg ) ((b + 2 c)g — xg') h 3 = ($,10) 


where primes denote differentiation with respect to to the 
new variable x. Note that the partial differential equation 
© has been reduced to the ordinary differential equation 
m- Equation ([Toll is difficult to solve in general. To 
demonstrate an exact solution it is necessary to make 
assumptions on the parameters and arbitrary functions. 

To progress we make the assumption 


g(x) = kh{x). 

Then equation m reduces to the form 

2ckx 2 [ck — ax b+c ] hh" 

3 a 2 x 2(b+c) - Aackx b+C + c 2 k 2 


( 11 ) 
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where f(t) is an arbitrary function of t. These symme¬ 
tries generate invariants that can be used to reduce the 
partial differential equation © into ordinary differential 
equations for further analysis. Using the symmetries in 
© in turn, or taking any linear combination, may be 
helpful in reducing the master equation into ordinary dif¬ 
ferential equations. Since the symmetries do not have a 
nonzero Lie bracket relationship we do not consider the 
optimal system in this paper. As all combinations of 
these symmetries can be contained in a general linear 
combination, we take 


aG\ + bG 2 + cG 3 
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Equation m is a nonlinear second order differential 
equation. We can reduce (fl2l) to first order if we let 
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Then the transformation (fT51) enables us to write ( 1121 ) in 
the form 
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where a, b and c are arbitrary constants. This combina¬ 
tion gives the invariants 



B = g(x)r b / c , (9c) 

where a ^ 0 and c ^ 0 . The quantities g{x) and h(x) 
are arbitrary functions associated with the Lie symmetry 
generators and arise from integration. 


Observe that (fTTl) is a Riccati equation in the variable 
y. Riccati equations can be transformed to second order 
linear equations. We let 


u(x) = exp 
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Note that in equation (fT5l) the term is the 

coefficient of the quantity y 2 in (fbTl) . Using m equation 
m is transformed to 
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where 


Hence we obtain the potentials 
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Therefore we have the remarkable feature that the second 
order nonlinear equation (1121) has been transformed to 
the linear equation (fTfil) via the transformations (1131) and 

(USD- 

Particular choices of the constants a, 6 and c allow us to 
integrate ED- (ED or (fldl) and obtain potentials for the 
gravitational field. There are several classes that arise 
which we consider in turn. 


IV. CLASS I 


In this class we make the restrictions that a / 0 and 
c ^ 0. These conditions arise so that the invariants © 
exist. 
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where x = exp ^ J* —j ^ ''j r 1//c and n is an arbitrary con¬ 
stant of integration. We believe that this solution is not 
contained in the literature. 

It is possible to evaluate the integral containing x in 
(121 all when m = 0. In this case 
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and h is a new constant. For this special case the explicit 
dependence on the variable x is fully specified; there is 
freedom only in the variable t. 


A. Case 1(a): 6 = 0 


B. Case 1(b): b = -c 


If we set 6 = 0 then the coefficient of the linear term 
in u in ED disappears. In this case we have 

u" + [3a 2 a: 2c (2ax c — k) + c 2 k 2 (k — 6 ax c ) 
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which is a separable equation. This can be easily inte¬ 
grated to give 
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where m = ^ is a new constant. 


If we set 6 = — c then equation ED becomes 
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which is a simpler form. It is interesting to note that this 
case produces the Euler equation (1^1) . We can integrate 
(l23l) to obtain 
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Hence we have the potentials 
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where m = Cic 3 3cfc and n = C 2 C 3 2c * are constants. 
This is another new solution to the master equation. 


C. Subclass 1(c): b = —c, k = —3a/c 


Note that in (l26l) . k ^ —3 a/c. With the values b = 
—c and k = —3a/c the transformation (fl5l) leads to an 
inconsistency. This means that we have to integrate CUD 
or (TH1) for this case. If we set k = —3a/c and b = — c, 
then equation (fl2l) becomes 

24 x 2 hh" - 2Ax 2 h' 2 + 24a :hh' + c 2 h 2 = 0, (27) 


which is greatly simplified. Now m can be integrated 
to give 
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where m and n are constants of integration. 
Hence we get the metric functions 



This is also another new solution to the master equation 

©■ 


where the constants d and rh result from the integration. 
Then equation (fl5l) yields 
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where m = ^ is a constant. 

Hence we can generate the potentials 
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where x = exp ^ J* —c 1 ^ c and n is new arbitrary 
constant. This is a new model. 

This class of solution contains a previously found 
model. To show this we set 


a = c = k = 1 and f(t) = t. 
Then solution (l33l) becomes 
A = B 
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which is a self-similar solution for the master equation. It 
has been obtained previously by Abebe et al. [23| when 
analyzing a conformally flat radiating star. 


D. Case 1(d): b = —2c 


V. CLASS II 


If we set b = — 2c then the coefficient of the linear term 
in u in © disappears as in Case 1(a). Then we obtain 
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In this class we set a = 0. Then the invariants © do 
not exist. For this category of solution we then have to 
consider the symmetry 
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The invariants (1571) reduce the master equation © to 

2c 2 gg'h' — c 2 (2 gg" + g ,2 ) h + 2 c(b + c)g'h 2 

+ (3 b 2 + 6 be + 2c 2 ) h 3 = 0. (38) 

Here the primes stand for the derivatives with respect 
to the independent variable t. Equation (1551) is an Abel 
differential equation in h. It is not possible to integrate 
it in general. Particular solutions do exist as we now 
demonstrate. We set 

/ ,/q \ 

(39) 

to simplify the Abel equation. Then (1551) becomes 

6 gg'h! — 3 (2 gg” + g' 2 ) h ± 2\^3g'h 2 = 0. (40) 

The advantage of (1401) is that it is a Bernoulli equation in 
h. It is interesting to note that it can be integrated even 
when the function g is unspecified. We integrate (l40l) to 
obtain 


Using the invariants (l44l) the master equation ([6]) is re¬ 
duced to give 
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where the primes denote derivatives with respect to the 
independent variable r. Equation (1451) is highly nonlinear 
and cannot be easily integrated in general. We assume 
that 
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so that 651 becomes 
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We observe that (1571) may be treated as a quadratic equa¬ 


tion in the quantity 4-; this gives 
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The two roots in (1451) can be integrated to yield two dif¬ 
ferent classes of solutions 
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which is a new exact solution for the shear-free model. 

Other exact solutions to (1551) exist but they may not 
be realistic. For example, (1551) can be integrated to yield 
exact solutions if we assume c = —6 and g(t) = kh(t). 
Unfortunately the model then become unphysical as the 
heat flux vanishes when g(t) = kh(t) and both the 
tangential pressure and the energy density vanish when 
c = —b. 


VI. CLASS III 

In this class we set c = 0. Consequently the invariants 
© are not defined. Therefore in this case we consider 
the symmetry 

aGi + bG 2 = (b — + B-gg + a /W^> (43) 

from ©. Note that the invariants arising from (1451) are 
given by 
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where a(r) = 2ay/a 2 + abkr + b 2 k 2 r 2 and d is an arbi¬ 
trary constant of integration. 

Hence we have 
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which are two classes of solution to the master equation. 
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This category of solutions reduces to known models. 
We can show this by setting 


f(t) = 1 and b = d = k = 1. 
Then the solutions (15(1 and (1CT1) become 
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- 4/3 f 1 a ( r ) 
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(52) 
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respectively. These solutions were previously obtained 
by Abebe et al. (23| for a radiating star which has the 
property of conformal flatness. 
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where if = r x / c exp ck . The matter 
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— 6 amnip)\ ck 3 (a — ck)(3a + ck) (mip + n) 2 


(55a) 


P II 


P± = 


Ack 

x (ip -1 / 2 (mip + n)j 3a+ ° 

= [2a (Gamnip — ck (m 2 ip —4mmp + n 2 ))] 
k 2 (ck — a) (3a + ck) (mip + n ) 2 


x (ip 1//2 (mip + n)j 3a+ck , (55b) 

. _ 4 cfc 

(a + ck) (ip - 1 / 2 (mip + n)J a+ ° (12 ammp 


Q = 


+ck(mip + ny J ck 3 (3a + ck) (mip + n)‘ 

Ack ~\ 1 

x (ip - 1 / 2 (mip + n)j 

k 
r 


— (mip 1 / 2 + nip 
r V 


- 1 / 2 ) 3 ^ 


-1 


(55c) 

(55d) 


VII. EQUATIONS OF STATE 


From the above we generate the linear barotropic equa¬ 
tion of state 


In relativistic astrophysics it is important that the 
model should admit an equation of state on physical 
grounds. Many exact solutions for a radiating star that 
have been found before do not satisfy this condition. We 
can report that particular classes of models in this paper 
do admit a linear equation of state which are barotropic. 
We give each of these cases below and associate them 
with the relevant Lie symmetry. 


P\\ = A/q A = 


provided that k ^ — —. 


ack 


v 2 - ack + c 2 k 2 ’ 


(56) 


(a) The generator aGi — cGi + cGs (k A —^) 

Using the generator aGi — cGi + cGz we can obtain a 
model from (1261) which admits an equation of state. The 
line element for this case is 


(b) The generator aGi — cG 2 + CG 3 (k = — —) 


The generator aGi — CG 2 + cGz helps to find another 
model (EH1) when k = — —. This model also admits an 
equation of state. The line element for this case is given 
by 


ds 2 = 


(mip 1 / 2 + nip 1 / 2S j 


\ 2 ck 

2 

1 

j 3a+ck 

r 

L/(*)J 


dt 2 


+ 


r fc n2 


[dr 2 + r 2 (d 6 2 + sin 2 ddcp 2 )] j , (54) 


ds z = 


mp 


e x P (f| log 2 v) 


-(1 /f(t)fdt 2 


x [dr 2 + r 2 ( dd 2 + sin 2 dd(p 2 )]\ , 


■ 9 (a/(cr)) 2 
(57) 







































where p = r 1 / c exp 
come 


(/ af(t )) 


. The matter variables be- 


/x = [13 (c 4 log 2 tp — 48toc 2 logy? + 24 (c 2 + 24m 2 )) 

1 


x exp (^-c 2 log 2 <p 
P\\ = [(48mc 2 logy? — c 4 log 2 ip — 24 (c 2 + 24m 2 )) 
x exp (ic 2 log 2 y; 

p_L = [(48?nc 2 log p — c 4 log 2 p + 48c 2 + 576m 2 ) 

1 


[2592a 2 nV m ] , (58a) 


[864a 2 nV m ] , (58b) 


x exp ( —c log p 


[648a 2 nV m ] , (58c) 


9 = 


3a 


rup 


cr ex p (fl lo § 2 <p) 


-i 


p ii- 


(58d) 


We observe that this case also admits an equation of state 


P\\=Xp, A = — —, 


(59) 


which is linear and barotropic. 


(c) The generator 6 G 2 + CG 3 

The generator 6 G 2 + cG% can also associated with an 
equation of state. The line element for this model be¬ 
comes 


d S 2 = _^ g ’ 2 r ±2^ dt 2 +g 2 r ±^-2[ dr 2 

+r 2 (d6 2 + sin 2 ddcj) 2 )] . (60) 

from (PS1) . In the above we have set the arbitrary con¬ 
stant d = 0 without any loss of any generality. The 
matter variables become 


d = 

14r T v^ 

(61a) 

39 2 ’ 


T 2 

4r + \/3 

(61b) 

P\\ = 

39 2 ’ 


r T ~75 

(61c) 

V .L = 

9 ’ 

9 

9 = 

r ±Ai-ii -1 

gr 3 p ||. 

(61d) 


This solution also satisfies the barotropic equation of 
state 


P\\=Xn, a = --, 


(62) 


VIII. DISCUSSION 


Our treatment indicates an interesting interplay be¬ 
tween Lie symmetries, Lie algebras, nonlinear differen¬ 
tial equations and a radiating star in relativistic astro¬ 
physics. We have analyzed a radiating star with a shear- 
free matter distribution with anisotropic stress using the 
Lie analysis of differential equations. A systematic study 
of the master equation, governing the evolution of the ra¬ 
diating star, was undertaken. Three classes of new exact 
solutions were generated. In each case the Lie infinites¬ 
imal generators can be identified and the gravitational 
potentials may be written explicitly. In certain cases we 
find that our models contain solutions found previously. 
In particular we regain the conformally flat solutions of 
Abebe et al. [23|. We also find that three exact mod¬ 
els, for specific forms of the Lie point symmetry, contain 
a linear equation of state. The line elements and the 
equations of state can be written in explicit form. We 
summarize our results in Table Q] and Table El In Table 
[2 the Lie point symmetries and corresponding invariants 
are listed. The restrictions on the arbitrary functions 
g(x) and h(x) and the constants a, b, c, and k are given. 
The resulting gravitational potentials A(r,t ) and B(r,t ) 
are written explicitly. Table |TT] identifies the symmetries 
with those particular exact solutions for which an explicit 
linear barotropic equation of state exists. 

The solutions found in this paper may be used to study 
the physical properties of a relativistic radiating star. To 
illustrate this we consider the temperature profiles. The 
causal heat transport equation becomes 


T (qB) t + ABq = --^(AT) r . (63) 


The coefficient of thermal conductivity k and the relax¬ 
ation time r are chosen as 


K = jT 3 t c , 




(64) 

based on physical grounds as shown by Govender and 
coauthors [HI, H3]- In equation (RTH) a > 0,/3 > 0, and 
7 > 0 are constants and r c is the mean collision time 
between massless and massive particles. We consider the 
special case of mean collision time a = 0 for simplicity. 
Then (HHIl) can be integrated to give 


(AT ) 4 = -- 
a 


/3 j A 3 B(qB) t dr + J A 4 qB 2 dr^j + w(t ), 


(65) 

which is the causal temperature. If we choose the grav¬ 
itational potentials in (1421) . the causal temperature can 


which is linear. 














9 


written explicitly explicitly as 


rjn 4 


(w(t) — r^g'(t ) 4 ^30d/3 + y/3ar^g(t ) 3//2 
+24/3-\/g(tj'j 4 ag(t) (d + y/g{t)j 
16 (d + ^/g(J)^ 4 


9 r^g{t) 2 g’{t) A 


( 66 ) 


Figure ID then gives the graphical behavior of the temper¬ 
ature when a = d = —1, g(t) = t = 0.5 and w(t) = —1. 

We find that the temperature is decreasing as we 
approach the boundary and the causal temperature is 


greater than the Eckart temperature. This is consistent 
with other treatments (see for example the recent analy¬ 
sis of Reddy et al. 0). 
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FIG. 1: Causal temperature (dashed line), noncausal temper¬ 
ature (solid line) versus r 


TABLE I: Symmetry, invariants, restrictions and resulting gravitational potentials 


Symmetry Invariants 


Restrictions 


Gravitational potentials 


aGl+bG2 X =^0S1„ 

+cG 3 A = $$ 


a n f rx Sacx c 1 exp (ax c -ck) ~ . 

f(t) ^ P exp(^|j^-) (6a 2 :r 2c —10acfca; c +3c 2 fc 2 )+ra ^ ^ 

B = - k f( t '> a where x = exp ( f* —rrpr) r _1 ^ c 

_ ip(/‘Jy) , \ afW) _ 

( 8 ac exp ( 2ax , — ^ (ckx c — a) 

f X “7- /, _ c \ r V J -T-7 

x yexpy 2q, ^ fc —J [6a 2 +3c 2 fc 2 a; 2c — 10ackx c ]-\-mx 2c J 


f * cdt 


x ex P (f* Tf(t) 

I r* _6dt_\ 

J af(t) J ’ 

B = g(x)r b/c 


g = kh, b = 0 


t cdt 

fit) 


dx 


_ ft cdt \ 

J af (t) J 


B = 


k (p-cdt) A > where j = exp (f t ^i)) r 1/c 
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TABLE II: Equation of state 

Symmetry Gravitational potentials Equation of state 






